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Abstract 

The way particles interact with turbulent structures, particularly in regions of 
high vorticity and strain rate, has been investigated in simulations of homogeneous 
turbulence and in simple flows which have a periodic or persistent structure e.g. sep- 
arating flows and mixing layers. The influence on both settling under gravity and 
diffusion has been reported and the divergence (compressibility) of the underlying 
particle velocity field along a particle trajectory has been recognized as an impor- 
tant quantity in quantifying these features. This paper shows how these features 
can be incorporated in a formal way into a two-fluid model of the dispersed parti- 
cle phase. In particular the PDF equation for the particle velocity and position is 
formerly derived on the basis of a stochastic process that involves the statistics of 
both the particle velocity and local compressibility along particle trajectories. The 
PDF equation gives rise to contributions to both the drift and particle diffusion co- 
efficient that depend upon the correlation of these quantities with the local carrier 
flow velocity. 
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1. INTRODUCTION 

There are two aspects of the motion of particles in turbulent flows that have not been 
properly incorporated in a rational way into a two-fluid model of a dispersed particle 
flow, namely the influence of persistent structures in the underlying carrier flow, and the 
occurrence of drift (either under the influence of gravity or as a result of inhomogeneity 
in the underlying turbulence). In their numerical simulations of particle settling in homo- 
geneous turbulence and in cellular flow fields, Maxey and his co-workers have shown for 
instance that turbulence can enhance the settling of small particles, (Maxey & Corrsin 
1986, Maxey 1987, Wang & Maxey 1993). In particular Maxey (1987) showed that in 
situations of weak particle inertia (i.e. particle relaxation times ^ the typical time scale 
of the turbulent structures in the flow) the net settling velocity of an ensemble of 
particles in a homogeneous flow fleld was related to its value in quiescent flow by the 
relationship 



V = - 



f (u(x,t)V-Vp(Xp(x,t I s),s))rfs, (1) 

^ 



^Present address: School of Mechanical & Systems Engineering, Stephenson Building, Newcastle Uni- 
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where ( ) is an ensemble average; u(x, t) is the carrier flow turbulent velocity field at 

position X for times t integral time scale of the turbulent motion; V • Vp (Xp(x, t | s), s) 
is the divergence of the particle velocity field Vp (y, s) with respect to the spatial position 
y measured at y = Xj,(x, i | s) at time s where Xp(x, t\s) is the position of a particle at s 
which arrives at x at time t. The particle velocity field Vp (y, t) is defined as the particle 
velocity field arising from one realisation of the fiow field u (y, t) with a prescribed set of 
initial conditions at s = for the particles which are the same in each realisation of the 
flow fleld. 

The divergence of the particle velocity field is a measure of the local compressibility of 
the particle fiow. The presence of gravity means that the particles move in a preferential 
direction which in turn means that the correlation of the fiuid velocity with the 'local' 
divergence of the particle fiow field is non-zero. For the case when the particles almost 
followed the flow, Maxey was able to relate the local compressibility of the particle flow 
fleld to the local straining of the underlying carrier flow and showed that the value of 
the correlation would lead to an enhancement of the gravitational settling. Subsequently 
Wang & Maxey(1987) explained this result in more detail by looking at the way particles 
move around the edges of vortices; in particular their results could be explained by the 
streaming of particles between vortices which always lead to an accumulation of particles 
on the down flow side of vortices (i.e. in the direction of gravity). They referred to 
this process as preferential sweeping. This however is not a unique result. For instance 
depending upon the particle Froude number, Davila and Hunt (1999) have shown it is 
possible for the opposite to occur. 

Whatever the particular route the particles take through a flow field (with or without 
gravity), the compressibility of the particle flow fleld measured along a particle trajectory 
is an important consideration in the way we assess the influence of structures. The 
work presented here shows that Maxey's expression for the drift is a much more general 
result appropriate for inhomogeneous as well as homogeneous fiows with or without the 
presence of gravity. Indeed it is shown that the compressibility of the particle fiow can 
influence not only the drift but also the particle dispersion. As a prelude to the full 
two- fluid formulation I first consider in Section 2 the analysis of particle dispersion and 
drift in a compressible fiow field in which the statistics of both the particle velocity and 
the divergence of the particle fiow field along a particle trajectory are prescribed and 
correlated. Then finally in Section 3 these features are incorporated into a two-fiuid 
model of the dispersed particle phase, based on the so-called pdf approach - the focus 
here being on the derivation of the appropriate transport equation for the particle phase 
space probability that a particle has a velocity v and position x at time t. Some of the 
features of this analysis are illustrated in particle dispersion in a random array of counter- 
rotating vortices in both homogeneous and inhomogeneous situations depending upon the 
prescribed statistics. 

2. PASSIVE SCALAR DISPERSION IN A COMPRESSIBLE FLOW 

2.1 Gaussian and non-Gaussian Lagrangian Statistics 

Analyses of this sort have been done before for passive scalar diffusion in incompressible 
flow in which case only the statistics of the particle velocity along a particle trajectory 
are required. In the case of a compressible flow, moments associated with the process 
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[vp (s) , V ■ Vp (s) , s G t] appear as a natural consequence of the transport and the com- 
pressibihty of the flow, where both the particle velocity and the divergence of the particle 
velocity fields are measured along a particle trajectory. 

The starting point of the analysis is the conservation equation for the particle mass density 
p(x, t) at position x = [xi,a;2,a;3] at time t, namely 



dp 

di 
Dp 



-V-{vpp(x,t)} 



or — = -pV-Vp(x,t). (2) 

where is D/Dt rate of change along a particle trajectory. Given some initial distribution 
p(x, to) at time t = 0, the solution is formerly 

p(x, t)=p (Xp(x, t|0), 0) exp |- ^* V ■ Vp (Xp(x, t\s), s) rfsj . (3) 

where Xp(x, t\s) is the position at time s of a particle arriving at x at time t. Given that in 
principle we can define a particle velocity field Vp(x, t) for any realization of the underlying 
carrier fiow field, then the problem of particle dispersion and settling is identical to the 
problem of passive scalar dispersion in a velocity field differing only from the normal 
case considered in that the particle velocity field is compressible rather than solenoidal. 
Replacing Xp(x, t|0) by x — Jq Vp (Xp(x, s) ds in Eq.Q we obtain 

p(x, t) = p ^x — y Vp(x, t I s) ds, 0^ exp ^~ J V ■ Vp(x, t \ s) (isj . (4) 

where Vp(x, t \ s) and V ■ Vp(x, t | s) are used as shorthand notation for the explicit values 
of the particle velocity and divergence along particle trajectories that pass through (x, 
namely 

Vp(x,t I s) = Vp(Xp(x,t|s),s) V- Vp(x,t I s) = V- Vp(y,s)]y^Xp(^_^l^) (5) 

We shall sometimes abbreviate these quantities still further to Vp(s) and [V ■ Vp] (s) respec- 
tively. By making certain assumptions about the statistics of the process [vp(s), V ■ Vp (s)] 
for < s < t , then we are avoiding the non-linearity of the diffusion process that is im- 
plicit in the relationship between Lagrangian and Eulerian timescales. As a result it is 
shown in Appendix A that if this process is jointly Gaussian then the particle drift ve- 
locity is given precisely by the term in Eq([T]) and the diffusion coefficient consistent with 
Taylor's theory. Explicitly the particle mass current is given by: 

(pvp(x,t)) = |(vp(x,t))-y*(v;(x,t) V-Vp(x,t |s))rfs|(p(x,t)) 

- ds (v; (x, t) v; (x, t I s)) ■ V (p(x, t)) (6) 

where (x, t) is the ffuctuating part of Vp (x, t) relative to its mean. The first bracketed 
term on the RHS (the drift term) in this equation is identical to the drift term derived by 



^It is implicit here that the divergence be appHed to the spatial components of the particle velocity 
field and is not meant to operate on x in the vector function Vp (x, t \ s) . 
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Maxey if we substitute V°for (vp(x, t)). However with the assumption that the statistics 
for the underlying particle velocity field are Gaussian, we end up with a more general 
result which includes a gradient diffusion flux. If in general the statistics of the process 
[vp(t), V • Vp (t)] are non-Gaussian then it is shown in Appendix A to first order in the 
triple moments of the process, that the particle mass current is compounded of a drift 
term 

Vd = (vp(x, t)) - dsi (V ■ Vp(x, t I si)Vp(x, t)J> + 

+^ (v;(x, t)W ■ Vp(x, t I Si) V • Vj,(x, t I S2)) (7) 

and a gradient diffusion term with diffusion coefficients Dij 

Dij = J dsi (?/p,(x,t I si)Vp.(x,t)) - J dsi j ds2 (up. (x, t)t/p.(x, sajV • Vp(x,t | s^)) 

(8) 

2.2 Comments on the process [vp(s), V • Vp ; < s < 

It is important to recognize that the statistics of the process [vp(s), V ■ Vp(s) ; < s < 
which we will call [q(s)] for short, does not depend on the initial concentration. If it 
did, then its statistics would be related to the particle density weighted averages we are 
trying to calculate in the first place . How these statistics are obtained is clear: a particle 
trajectory is soved backwards in time starting from x at time t using the values of the 
particle velocity Vp(s) along its trajectory which in turn are derived from the statistics of 
the velocity field Vp(x, s) where < s < 1 No restrictions are placed on the point the 
trajectory goes through at time zero. In the actual problem of interest we might want to 
know the average particle velocity at x at time t knowing say that the particles started 
out at xq at time zero. So these particular particles will choose a particular subset of the 
statistics of the process [q(s)] in arriving at x at time t. That is, we are selecting only 
those trajectories of all those trajectories defined by the process [q(s)] that go through Xq 
at time zero from which we could compute the particles average velocity at x at time t. 
Put another way, we are trying to evaluate the particle statistics from a set of statistics 
which are independent of where the particles start from in the actual problem of interest. 
You can see this more transparently in the way the concentration is calculated. You start 
off with some prescribed statistics for the process [q(s)] found by starting a test particle 
off at X and solving the equation of motion backwards in time from t. That is you solve 

= Vp(Xp, s) withXp(t) = X (9) 

backwards in time to find the values of Xp(0)and the value of exponential of the integral 
of the value V • Vp along a trajectory (which gives the fractional change in an elemental 
volumee at time t along the trajectory relative to its initial value i.e. the value of the 
elemental volume deformation J{t) = ^"^^^'g^^'* . You then calculate the concentration 
that particles would have at [x, t] if they started out at time with some concentration 
p(Xp(0),0) by multiplying this concentration by J{t). If the concentration at Xp(0) hap- 
pens to be zero, then the concentration at x is zero. The fact that there may not be any 
particles at Xp(0) doesn't affect the statistics of the process [q(s)]. The process just tells 
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you where you might find some particles at time but if there aren't any, then that's 
because of the initial conditions. The process [q(s)] doesn't know about initial conditions 
or concentration. It's entirely determined from the statistics of Vp(x, t) derived from some 
test particle at x at time t in the way we have prescribed. 

2.2 Dispersion in homogeneous staionary turbulence and comparison with Taylor's Theory 



It is revealing to compare these results derived for passive scalar diffusion in a compress- 
ible flow fleld with G I Taylor's classic theory for diffusion by continuous movements in 
a homogneous stationary flow fleld? We recall that Taylor's results are based on the as- 
sumption that in the limit of the dispersion time t — )■ oo this time can be divided into a 
large number of time steps (each step >> the integral timescale of the turbulence) so that 
the distance travelled in one time step will be uncorrelated with the distance travelled in 
the next. This leads to Gaussian statistics for the particles displacement. In particular 
the diffusion coefficient is written as 



oo 



A, = Ri^ds, (10) 

where Rij{s) defines the velocity autocorrelation <^Vp^(x, | Si)vp^. (x, | ^2)^ for which 
s = si — S2, and Vp-(x, | Si) and Vp^(x, | S2) as before are the velocities of a particle 
measured at time Si and S2 that particle starting out at some arbitrary position x at 
some arbitrary time t = 0. The important requirement is the distances at which these 
measurements take place are on average very far away from the point of release x so that 
the process [vp(x, | s)] is stationary i.e for this to occur s/T^ ^ 1 where is the 
Lagrangian integral timescale. The rate of change of the mean square displacement is 
given by 



2{Vp,(x,0 I «)X„(x,0|()) =2 /"{Vk(x,0 I «)vp,(x,0 | s)) ds (li; 





and with the assumption that one can replace the lower limit by say r such that t 
^ t — r ^ Tl, so that during the interval t < s < t,Vp.(x, | s) is stationary, one 
arrives at the classic result that 



|(^P.(x,0|t)^ 



- t— >oo 



Rij{s)ds (12) 



a result which is self consistent with a Gaussian or gradient diffusion process with a 
diffusion coefficient given by the RHS of Eq (fT2|) . Returning to the form for the diffusion 
coefficient defined in Eq.() that is the extension to non-Gaussian fields, we note that with 
the Taylor Gaussian assumption for t — )■ cx3, we are left with the result that 

A,=^*rf3i(v;(x,t |s)v;^ (x,t)) (13) 



which since <^Vp^(x, t | s)Vp^. (x, t)y is dependent only on t — s we get the same result as in 
Eq.([l2]) if 

'4(x,t 1 3)v;^(x,t)) = (v;^(x,o 1 3)v;^(x,o 1 1)) (m) 
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If the flow is homogeneous and stationary then these correlations since they refreed to the 
same particle measured at two different times, will be independent of labeling position 
and times. That is we could change the labeling time from t to in the correlation on the 
RHS of Eq.( fT^ and retain the labeling position x without changing the result. So the 
relationship in Eg. (1141) is valid. The relevance of this independednce on labelling times 
and positions is even more revealing when we consider the general result for the rate of 
mean square displacement given in Eq. dlll) for all t and compare it with the form derived 
from the continuity equation 

dp d 

on the form for Dij in Eq.([S]) appropriate for non-Gaussian fields. That is if we release 
particles at time and measure the dispersion at time t then the form of Dij would imply 
that 

x,0|i),t]Xp^(x,0|t)) = (vp^(x,t)Xp,(x,t|0))-|^*dsi f^ds2 (V-Vp(si)v;^(s2)v;^,(t))+. 

(15) 

The statistics associated with the correlation on the LHS of this equation is different 
from that determining the first term on the RHS of the equation: in LHS case, we have 
statistics derived from two Lagrangian variables where arguments for stationarity can 
only be invoked when t — t- oo, whilst the case of the RHS is derived from a Lagrangian 
and an Eulerian variable. The two are only equal when t — )■ oo or at small times t <^ 
when the second term on the RHS is 0{t/Ti) smaller. The term on the RHS is clearly a 
measure of the difference in the two sorts of statistics. 



3. PDF FORMULATION 



This represents an extension of previous work by this author (Reeks 1991, 1992) and 
several others (Zaichik 1991, Swailes 1997, Hyland et al. 1999, Pozorski & Minier 1999 
and Simonin et al. 1999) in using an equation for the particle phase space probability to 
formally derive the two-fluid continuum equations for the particle phase. 

2. 1 Definition and prescription of particle velocity field and its divergence 

Using Stokes drag as an example, the particle equation of motion can be written as 

d\ , , , (ix , , 

^ = ,fl{u(x.i)-v}+g ; -=v (16) 

where as before u(xp, t) is the underlying carrier flow velocity at position x at time t. The 
solution can be written in several ways. First solving the set as a time problem, 

v = Vp(u;,y,t/|t) ; x = Xp(u;, y, t' | t) (17) 

i..e. the solution is the particle velocity/position at time t, for a particle with initial 
velocity uj and position y at time t\ allowing for the possibility of t' being in the past or 
the future in relation to t. Clearly these functions deflne the inverse relations 

a; = Vp(v,x,t I t') ; y = Xp(v, x, t | t') (18) 

So using these equations we could eliminate y from Eq. f|T7j) and write an alternative 
solution, namely 

V = Vp(a;,x,t' I t). (19) 
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Figure 1: Pairs of counter-rotating vortices generated from random symmetric shear flow 



That is the particle velocity at position x at time t given that the particle velocity started 
out at time t'with a velocity u. \p{(jj,^,t' \ t) is the particle velocity field (in the context 
of the passive scalar dispersion in Section 2), which satisfies the equation 



Vp(a;,Xp,t'|t). 



So we have 



dt 'P^-'^-P' 
d \ dXpX<^,y,t'\t)\ ( dvp 



dt 



dX. 



Given that the initial conditions imply that 



dXp. 



at t — t ^ 



then we would have directly from Eg. (1211) that 
9Xp(u;,y,t'|t) 



'dXp 



J 



dy 



exp ds V ■ Vp(c<;, y, t'ls) 



ly=Xp(x,t|s) 



(20) 
(21) 

(22) 



(23) 



2.2 Closure of the PDF Equation 

If H^(x, v,t) is the phase space density for a particle with velocity v and position x at 
time t subject to the equation of motion defined in Eq. fll6p for one realisation of the 
carrier fiow filed u(x, t), then the equation for (Vr(v, x, t)) the PDF for a particle to have 
(v, X, t) is obtained by averaging the Liouville equation thus, 



(M') = -|;-/'{u'(x,i)lV) (24) 



where (u(x, t)) and u'(x, t) are the mean and fiuctuating components of u(x, t). We 
require therefore a closed expression for (u'(x, t)Vr). In reality we consider a closed 
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expression for the specific case when is a response function G, that is it is the solution 
for an instantaneous point source 5(v — v')5(x — x.')S{t — t'). Thus (G) is the solution of 
the PDF equation Eg. (I24|) with the instantaneous point source added to the RHS of the 
equation. Knowing (G) we have for (iy(v,x, t)) 

(ly (v, X, t)) = J dv rfx(G(v, X, t|v', x', t'))p(v', x', t')dVdx.' (25) 

where p(v', x', t') is some initial distribution of (W) at time t'. With reference to Eqs. fll7l) . 
we can write the solution for (G) formally as 

(G) = (5(v - Vp(v', x', t' I t))5(x - Xp(v', x', t' \ t)) (26) 

However using the definition of the particle velocity field Vp(v', t'|x, t) we can write this 
alternatively as 



(G) = {6{v - vp(v', t' I X, t))5(x - r ds vp(v', t' I y, s) 

Jt' 



exp{~ dsV ■ Vp(v', t' I y, s] 



}) (27) 

y=X(x,t|s) 



Similarly we can write down formally an expression for (Gu). This expression together 
with that for (G) are in form that we can process in a similar manner to the evaluation of 
(pvp) for the passive scalar case: the only difference here is we are considering a process 
[vp(s), V ■ Vp(s), u(s) ; t' < s < t] as opposed to [vp(s), V • Vp(s) ; t' < s < t\. We show in 
Appendix B that if this process is Gaussian, then f3{Gu) is given exactly by 

^(u(x, t)G(v, X, t| v', x', t') = - (^/x • ^ + A ■ ^ (G) + 7(G) (28) 

where 

fi = /3(u'(x,t)vp(t)) 
A = /3(u'(x,t)xp(t)) 

7 = -/3 fds (u'(x,t)V-Vp(s)), (29) 

where 

Vp(s) = Vp(v,t'|Xp(x,t|s),s) V-Vp(s) = V - Vp(v',t'|y,s)ly=xp(x,t|s) (30) 

Xp(t) = Xp(x,t|0) = /* Vp(v',t'|Xp(x,t|s),s)rfs (31) 

Jo 

The form of the net force per unit mass of particles due to the turbulence given in 
Eq. (l28|) is therefore composed of two parts: a diffusive force (gradient of a stress tensor) 
which depend upon gradients in the particle velocity and position [the bracketed term 
on the RHS of Eq. (l28p ] and a body force which depends upon the local compressibility 
of instantaneous particle velocity field along a particle trajectory [the second term in 
Eq. (l28|) ]. The general form of this turbulent force has been obtained before by several 
authors (Reeks 1992, Swailes et. al. 1997, Pozorski and Minier 1999, Hyland et. al. 1999) 
but the precise form for the body force is different from the one derived here and leads to 
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Figure 2: Particle concentration profiles in random pairs of counter-rotating vortices 
(see FigHD 

the so-called problem of spurious drift; that is there is a drift term that persists in cases 
where the particles follow the underlying incompressible flow 0) so that where 

at equilibrium the particles ought to be fully mixed with the flow, the existence of the 
spurious drift leads to a build up of concentration in regions of low turbulence intensity. It 
is a feature that is common in certain types of simple random walk simulation of particle 
dispersion in inhomogeneous turbulence (where the underlying flow filed is essentially 1-D 
and cannot of its own accord satisfy continuity of flow if it is spatially varying. The form 
derived here does not suffer from this serious defect, the drift velocity in this case 
clearly vanishes when the particle follow the flow, because V ■ Vp(s) is the same as 
that of the underlying carrier flow which is necessarily zero. 

As an illustration of the influence of turbulent structures I have considered the dispersion 
of particles in a random flow field which consists of pairs of counter-rotating vortices (see 
Fig.l) with randomly generated vorticity that shifts randomly in position as the timescale 
of the vorticity changes randomly from one value to the next in time. In the case of a 
flow field in which the location and periodicity of the structures is fixed, particles ac- 
cumulate at the stagnation points. That is the process is equivalent to diffusion plus a 
drift directed towards the stagnation point. As an example Fig.2 shows the difference in 
behaviour between a 1 D flow field in which the particles are constrained to move only in 
the x-direction and when they allowed to move in the y -direction (fully 2D vortex flow 
field). The difference illustrates the difference between spurious drift (arising from a 1 D 
carrier flow field which is compressible) and the case of an incompressible 2D carrier flow. 
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APPENDIX A 



Al. Gaussian Process 

We expand p (^r— /q Vp(s)c/s, 0^ in Eq.(j4]) as a Taylor's series about p(r, 0) so that formally 



p(r, t) = exp ■ 



/ V-\'p(s)ds+ / dsvp(s) 
Jo Jo 



d_ 
dr 



p(r,0). 



(32) 



We now suppose that both Vp(s) and V-Vp(s) to be a continuous processes whose statistics 
are correlated. That is Vp(s) is the limit of the discrete process 



=N~,oo [Vp(si), Vp(si)..Vp(Sj).., Vp(sAr)] 

jr with Nt = s. 



Similarly for V ■ Vp(s). For convenience we specify a vector q(s) 

q(s) = [Vpl (s) , Vp2 (s) , Vp3 (s) , V ■ Vp(s)] 



(33) 



(34) 



whose statistics we specify through the characteristic functional M [4>{s)] given formally 
by 



M[0(s)] 



exp I 



(s) ■ q(s)(is 



where 0(s) is an arbitrary vector function of time 



(35) 
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and we further assume that q(s) is Gaussian so that 
M [(j){s)] = exp (i f (q(s)) ■ (f){s)ds 



~\lo L ^^2(9i(^l)9i(^2))0i(Sl)0j(s2) 

where 

{{q^{s)) = ■ qiis,) = q, (si) - W.} 2 < 3 ; {{q,{s)) = ; q[{s) = q,{s)} i = 4 

We recognize from the definition of the characteristic functional that 

(p(r,t)) = M[^0(t)]p(r,O) 
d 

with (piit) = — — fori < 3 
dxi 

= 1 i = A 



(36) 



and 



(v.(t)p(M))^-^^^|^p(r,0)M.<3 



(37) 



(38) 



Substituting the Gaussian functional for M given in Eq. (l36|) into Eq (l38|) and performing 
the functional differentiation we obtain 



(vpfe(i)p(r,t)) = |(gfc(t)) + z {qA{s)qkit)) Ms)ds + ds {qi{s)qk{t)) (pii-s) j {p{r,t)) 

(39) 

Substituting the values for (j)i{s) in Eq. (l371) into Eq. (l39|) we obtain finally 
K(t)p(r,t)) = {(vp,(r,t)) - (vp,(t) V ■ v,,is)) ds^ (p(r,t)) 



-/;^^(vkW<(-))^(p(r,t)). 



dxi 



(40) 



Al. Non- Gaussian Process 

The same analysis can be extended to consider dispersion and drift in which Vp(s) 
and V ■ Vp(s) are jointly non-Gaussian in which we express the characteristic functional 
M [(f){t)] in terms of the cumulants of q (t) ,i.e. 



M [0(s)] = exp 



I 



\ 



(q(s)) -0(5)^5+ 



m=2 



^ /o '^■^1 /o ds2 /o (iSmX 



Ifji \ 771/ 





\ 


X 




\ 
/ 





(41) 



<ihisiWi2 (^2) --gL 



where 

M [0(s)] and Eq.i^ we obtain: 



represent the cumulants of q(t). Using this form for 



■m=l 



Q.'iSsi)..q[^ {Sra) q'k it) 
0ii(Sl)..0i^ {Sm) 



X 



(p(r,t)) (42) 
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with (j){t) given by Eqs( l37l) . So picking out the contribution to the drift and to the 
gradient diffusion from the correlation of the process [vp(s)] with the process [V ■ Vp(s)] 
we can write Eqs (H2]) more transparently as 



W, 



E 

m=l 



m=2 



I ||(v-vp(si)..v-vp(sjvp;(t)) 

■ Iodsi{vp/{si)vp/{t)) + 

l(..vp/(si)..v-vp(sjvp;(t))|| 




(p(r,t)) (43) 
9(p(r,t)) 



+ ...(44) 



So to first order in the triple moments of [q(t)] ,the convective velocity and Diffusion 
coefficients Dij are respectively 

Vrf = (vp(r,t)) - / dsi (V ■ Vp(si)vp(t)) + - / dsi rfss (V ■ Vj,(si) V ■ Vp(s2)vp'(t)) + ... 

Jo I JO Jo 

(45) 

Dij = / dsi{vpi'{si)vp/{t)) - / dsi / (is2 (V ■ Vp(si)..Vp/(s2)vp/(i(:)) (46) 
JO Jo Jo 

So the diffusion coefficient is derived from two parts: one which is appropriate for incom- 
pressible flows if the process and the other which is appropriate for compressible flows for 
non-Gaussian processes. 

APPENDIX B Evaluation of (Gu(x,t)) 
We can formally write Eq. (1271) as 



(G) 



exp 



d 



G(°)(v,x,t|v',x',t') 



(47) 

where v„(s) is used as shorthand for Vp(v', t/|X„(x, s) and a similar short hand of 



V ■ v„(s) for V ■ v„(v', t/|Xp(x, t|s), s). and \'„(s) is the fluctuating value of v„fs) with 



respect to its average value (vp(s)) . G'^'^''(v, x, t) is the response function 

= 5(v - (vp(t))5(x - f ds ((vp(.)) - x')) 

Jt' 

So as for the passive scalar case we consider the statistical process 

q(s) = (s) , V'p2 (s) , V'p3 (s) , U[{s), ^2(5), ^3(3), V ■ Vp(s), 



(4^ 



(49) 



with a given characteristic functional M[0(s)] which we will assume is a Gaussian func- 
tional. We have thus as before 



(G) = M[i0(s)] 
with (f)i{s) = 6{s - t)— + — for 1 < i < 3 

OVi OXi 




1 



for A < i < 6 
for 2 = 7 



(50) 



(51) 



12 



Particle drift in turbulent flows 



M. W. Reeks 



(„,(<)G) = f^£M£MG(.)(v,x,t|v'.x',t) (52) 

Performing this functional differentiation on the Gaussian Characteristic functional, leads 
to the closed expression 

(53) 
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